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I.  INTRODUCTION 


In  this  paper  we  derive  formulas  for  calculating  the  response  of  a  hollow 
cylinder  to  internal  loads  moving  with  uniform  velocity.  This  analysis  is 
based  on  the  historic  Green-Lame  theory  of  wave  motion1 2  and  recent  work  on 
free  vibrations  of  a  hollow  cylinder/  Only  steady  state  conditions  are 
considered;  transients  which  may  arise  at  the  beginning  of  motion  are 
ignored.  Formulas  for  moving  sinusoidal  loads  are  discussed  in  detail.  The 
response  to  a  moving  step  function  can  be  obtained  by  using  Cauchy's 
discontinuous  factor. 

The  transition  from  traveling  to  stationary  loads  is  obtained  by  a 
limiting  process.  The  biharraonic  functions  required  for  the  stationary  load 
arise  from  the  confluence  of  the  vector  and  scalar  wave  functions  at  zero 
velocity.  This  solution  is  also  derived  independently  from  appropriate 
elastic  potentials.  This  analysis  shows  that  the  calculated  displacements  for 
moving  loads  approach  the  calculated  values  for  stationary  loads  as  the 
velocity  of  travel  approaches  zero.  Physically,  this  requirement  is 
obvious.  Hence,  computations  for  low  velocities  will  serve  as  a  critical  test 
of  the  analysis  and  computations. 

The  results  are  required  to  interpret  strain  measurements  obtained  with 
an  instrumented  gun  tube  in  which  arrays  of  strain  gages  and  pressure  gages 
were  mounted  in  close  proximity  along  the  length  of  the  tube.  Tube  strains 
based  on  the  measured  pressures  were  calculated  according  to  thick  walled 
cylinder  theory,  using  appropriate  mechanical  properties  for  the  gun  steel. 

At  low  velocities  the  measured  and  calculated  strains  agreed  reasonably  well, 
but  at  high  projectile  velocities  the  calculated  strain  distribution  profiles 
were  significantly  different  from  the  measured  strain  histories. 

This  effect  is  well  known  in  the  theory  of  moving  loads/  Prof.  Ian 
Sneddon  of  Glasgow  University,  a  pioneer  in  the  field,  suggested  this  problem 
to  the  author  during  a  recent  visit. 

The  analysis  is  given  in  considerable  detail  to  facilitate  verification 
of  formulas  required  for  programming.  Formulation  is  limited  to  subsonic 
velocities  which  occur  in  practice;  analysis  of  the  moving  load  problem  for 
supersonic  velocities  is  quite  difficult  and  will  not  be  considered  at  the 
present  time. 
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1M.E.  Gurtin ,  Elasticity ,  Enc.  of  Phys. ,  Volume  VI  a/2,  Mechanics  II, 

Springer  Verlag,  New  York,  1972.  See  pages  212-214. 

2 

Armenakas,  G.  Herrmann  and  D.C .  Gazis,  Free  Vibrations  of  Circular 
Cylindrical  Shells^  Pergamon  Press,  New  York,  1969~. 

A.S.  Elder  and  K.L.  Zimmerman,  "Stresses  in  a  Gun  Tube  Produced  by  Internal 
Pressure  and  Shear,"  BEL  Memorandum  Report  No.  2495 ,  June  1975  (AD  A012765). 

L.  Taylor ,  Vibrations  of  Solids  and  Structures  under  Moving  Loads , 

Noordhoff  International  Publishing,  Gronigen,The  Netherlands,  1972.  See 
Chapter  17  and  Reference  203. 
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II.  FORMULAS  FOR  WAVE  MOTION  IN  CYLINDRICAL  COORDINATES 


In  the  absence  of  body  forces  the  equation  of  motion  may  be  written  in  the 
form* 


2* 

(A+2jj) VV»u  -  ji  V  x  (Vxu)  =  p  (1) 

3t 


* 

where  X  and  p  are  the  Lame  constants.  The  speed  of  the  dilatational  wave  is 

9 


Cj2  =  (A+2y)/p 


(2) 


and  the  velocity  of  the  shear  wave  is 


c 


2 


2 


=  w/p  , 


(3) 


where  p  is  the  density  of  the  elastic  material.  The  equation  of  motion 
becomes 


2* 

2  „„  2  „  \  3  u 

c.  V7*u  -  c„  V  x(Vxu)  =  — 

3t 


(4) 


The  Green-Lame  solution  is 


4 


u  =  V<j>-bVx\J; 


(5) 


A 

in  which  the  potentials  <f>  and  ^  satisfy  the  equations 


A  A. 

*the  notation  VV*u  is  to  be  interpreted  as  V(V*u). 

0 


L 


(6) 


2  u2a 

c,  v  <j> 


and 


(7) 


In  cylindrical  coordinates  we  have 


+  1  li 

r  3r 


+ 


2 

r 


(8) 


The  vector  potential  is  given  in  terras  of  three  components 


+  =  er  *r  +  e0^0  +  ez  i|»z  ,  (9) 


A  A 


where  er>  eQ ,  and  e are  unit  vectors  in  cylindrical  coordinates? 


l  A  A 


Z  2 

*  *  =  er  -  -J 

r 


1  a 
r 


2  36 


+  e0  v  *e 


[0  2  3 


2  36  *r 


(10) 


+  e  V  r 
z  z 


^  — - - - 

P.C.  Cfton  and  N.J.  Pagano ,  Elasticity:  Tensor ,  Dyadic ,  and  Engineering 
Approaches ,  £>.  Fan  Nostrand  Company,  Inc73  Prvnceton,  1967,  pages  245-265. 
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The  following  scalar  partial  differential  equations  are  finally  obtained: 


2„2  89 

c,  v  ♦  =  — J 

9t 


2„2  ^r  2  3 

C2  V  ^r  2  +  ’  2  30  ^0 

r  r 


3t 


2  rr 


2„2. 
C2  V  *0 


2_  3_ 

2  2  30  ^r 

r  r 


3t 


2  0 


2  u2, 

C2  V  ♦z 


3t 


2  Tz 


(ID 

(12) 

(13) 

(14) 


Four  components  of  the  displacement  vector  are  obtained  from  solutions  of 

these  equations.  Since  only  three  components  are  required,  we  can  specify  one 

additional  condition.  We  choose  a  relation  between  it  and  ill.  which  will 

Tr  r0 

facilitate  separation  of  variables,  following  the  ideas  outlined  in  Reference 

2. 

III.  SEPARATION  OF  VARIABLES  FOR  SINUSOIDAL  LOADS  MOVING  WITH  CONSTANT 
VELOCITY 


The  scalar  potential  can  be  obtained  directly  by  separation  of 
variables.  For  brevity  let 


Hj  =  cos(n0)sin[s(z-ct)] 


(15) 


H„  =  cos (n0)cos [s (z-ct ) ] 


(16) 


=  sin(n0)sin[s(z-ct)] 


(17) 
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=  sin(n9)cos [s(z-ct)] 


(18) 


and 


♦±  =  H.  f (r),  i  =  1,2, 3, 4. 


(19) 


In  these  equations  n  is  an  integer  since  the  stresses  must  be  periodic  in  a 
complete  cylinder.  The  constant  c  is  the  velocity  of  travel.  We  assume 


0<e<c<c2~e 


(20) 


where 


0  <  e  <  <  1 


(21) 


Let 


2  2  2 

«1  =  c  /cl  (22) 

8:2  =  1  -  «j2  (23) 


Then 


f"  +  ^f'  -  (\+  BjV)  f  =  0 
r 


(24) 
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and 


f  =  A  I  (0  sr)  +  A  K  (0  sr) 
ini  z  n  1 


(25) 


The  axial  component  of  the  vector  potential  can  be  treated  in  a  similar 
manner.  We  let 


*z  =  Hi8z(r),  i  =  1,2, 3, 4 


2  2.2 

a2  =  C2  'C 


2  2 
*2  "  1  "  “2 


(26) 


(27) 


(28) 


then 


g 


z 


=  A3In(02sr)  +  A4Kn(02sr) 


(29) 


The  radial  and  tangential  components  of  the  vector  potential  are 
coupled.  The  following  four  combinations  lead  to  separable  solutions: 


^r  ~  Hl8r *  ^0  H380 
^r  =  H28r *  ^0  H480 
^r  =  H38r *  *0  =  Hl80 


(30) 


(31) 

(32) 


*r  =  Vr’ 


^0  “  H280 


(33) 


The  combination  given  in  Eq.  (32)  is  discussed  in  Reference  2;  the  remaining 
combinations  are  introduced  to  satisfy  a  variety  of  boundary  conditions  on  the 
inner  radius.  We  note  that 


s 


3_ 

36 


H, 


-nH3 


(34) 


(35) 


3_ 

36 


H„ 


nH, 


(36) 


3_ 

36 


H4  -  nH2 


(37) 


On  referring  to  Eqs.  (12)  and  (13)  we  find 

r2gj!  +  rg'  -  (n2  +  1  +  &22s2r2)gr  +  2ngQ  =  0 
r2g*Q  +  rg'  -  (n2  +  1  +  022s2r2)gQ  -  2ngf  -  0 


(38) 


(39) 


If  we  let 


80 


"  8, 


(40) 
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we  find 


r2g^  +  rg'  -  [(n  +  l)2  +  e22s2r2]gr  =  0 


(41) 


and 


g 


r 


A5In+l(82sr)  +  A2Kn+l(S25t> 


(42) 


If  on  the  other  hand  we  let 


g 


(43) 


then 


r2gj!  +  rg'  -  [(n  -  l)2  +  322s2r2]gr  =  0 


(44) 


and 


Sr  '  Vn-l'V')  +  A8Kn-l'62sr) 


(45) 


In  this  report  the  radial  and  tangential  components  are  based  on  Eqs.  (40), 
(41)  and  (42). 

As  mentioned  previously,  a  complete  solution  to  a  specified  boundary  value 
problem  may  be  obtained  if  one  of  the  components  of  the  vector  potential  is 
set  equal  to  zero.  Hence,  although  unique  solutions  for  the  displacements  are 
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expected,  the  choice  of  potentials  leading  to  a  solution  is  not  unique. 

-  A  A  A 

Pao  and  Mow*5  use  vector  potentials  L,  M,  and  N,  which  are  derived  from 
scalar  functions  a,  \p,  and  x  in  their  analysis  of  diffraction  waves  in  an 
elastic  solid.  The  relation  between  this  solution  and  the  solution  of  Herrmann 
and  Gazis  may  be  obtained  by  equating  formulas  for  the  displacements. 

IV.  COMPONENTS  OF  THE  VECTOR  DISPLACEMENT 

Each  component  of  the  vector  displacement  consists  of  one  term  from  the 
scalar  potential  and  two  terms  from  the  vector  displacement,  according  to  Eq . 
(5).  These  terras  must  have  the  same  trigonometric  factor.  In  addition,  the 
signs  of  g^  and  g^  are  chosen  so  these  functions  involve  Bessel  functions  of 
order  n+l  when  gQ  =  -  g  .  The  sign  of  g  is  chosen  so  the  divergence  can  be 
written  in  the  form 


=  Hj^Cr),  i  =  1,2, 3, 4 


(46) 


where 


h(r) 


I  !_ 

r  3r 


(r*r) 


"  7  80  +  sg2 


(47) 


We  consider  four  cases,  as  shown  in  the  table  below. 


Pao  and  C.  Mow ,  Diffraction  of  Elastic  Waves  and  Dynamic  Stress 
Concentration  Crane  Russak}  Publishers,  T971 .  See  pages  217-239. 
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Case 


TABLE  1.  TRIGONOMETRIC  FACTORS  FOR  SCALAR  AND  VECTOR  POTENTIALS 


Function 

f(r) 

gr(r) 

8  0 

(r)  gz(r) 

h(r) 

I 

H1 

h4 

h2 

h3 

H4 

II 

h2 

h3  • 

~H4 

H3 

III 

H3 

h2 

"«4 

H1 

h2 

IV 

H4 

H1 

"H3 

"H2 

H1 

The  displacement  components  for  each  case  are  obtained  from  the 
potentials . 


Case  I 

ru  =  [rf'  +  srgft  +  ng  ]  H. 

b  z  1 


rv  =  [-nf  -  srgr  -  rg']  H3 

r  3 

rw  =  [srf  +  rg0  -  ngj  H2 


Case  II 


ru  =  [rf'  +  srg  -  ng  ]  H0 
b  z  l 


rv  =  [nf  +  srgr  +  rg']  H4 


(48) 


(49) 


(50) 


(51) 


(52) 
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(53) 


rw  =  t_  srf  +  _  rgQ 


ngr]  Hx 


Case  III 


ru  =  [rf'  +  srgQ  -  ngz]  H3 


(54) 


rv 


[nf  -  srgr  -  rg']  ^ 


(55) 


rw  = 


[srf  -  J7rgQ  +  ngr]  H4 


(56) 


Case  IV 


ru  =  [rf'  +  srgQ  +  ngj 


(57) 


rv  = 


[nf  +  srgr  +  rg']  H2 


(58) 


rw 


[-srf  +_  rg0  +  ngr]  H3 


(59) 


Next,  we  express  the  displacement  and  divergence  for  Case  I  in  terms  of 
Bessel  functions  of  order  n  and  n+1.  The  following  identities  are  used  to 
eliminate  derivatives : 


XI  '(X) 
n 


nIn(X)  +  XWX) 


(60) 


XK  '(X) 
n 


nKn(X)  - 


(61) 
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XI'n+1(X)  =  -  (n  +  l)In+1(X)  +  XIn(X) 

XK'  (X)  =  -  (n  +  1)K  ,,(X)  -  XK  (X) 
n+l  n+l  n 


(62) 


(63) 


On  carrying  out  details  of  the  analysis  we  find 


ru  =  {Ai  [nl^(3^ sr)  +  gjSrl^^sr)  J 


+  A2[nKn(S1sr)  -  gj srKn+1 (gjSr) ] 


+  A_nl  (g  sr)  +  A. nK  (g.sr) 
j  n  2  4  a  2 


-  A,.srln+i  (g2sr)  -  A^g^rK^^ (g2sr )  }  cos(n0)sin[s(z-ct)]  (64) 


rv  =  {-  A  nl  (g  sr)  -  A~nK  (g  sr) 

ini  2  n  1 


+  A3[-  nln(g2sr)  +  g2srln+1(g2sr)] 


+  A^[-  nKn(g2sr)  +  g2srKn+1 (g2sr) ] 


-  A^srl^j  (g2sr)  -  A^srK^j  (g2sr ) }  sin(n0)sin[s(z-ct)  ]  ,  (65) 


rw  =  {AjSrl^gjSr)  +  A^rK^gjSr) 
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A5srln(g2sr)  +  A^srK^B^r ) }  cos (n0) cos [ s ( z-ct ) ] 


(66) 


The  divergence  of  the  vector  potential  is 


V-*  “  [ (A3  +  A562)sl^(32sr)  +  (A^  -  A562)sKn(32sr) ] sin(nQ)cos [ s ( z-ct) ] 

(67) 


Formulas  for  the  displacement  and  divergence  corresponding  to  Cases  II,  III, 
and  IV  can  be  written  down  by  inspection* 

The  divergence  of  the  vector  potential  enters  into  calculations  of  the 
rotation  vector.  The  analysis  prior  to  Part  IX  of  this  report  follows 
Reference  2  in  using  a  vector  potential  with  nonsolenoidal  divergence,  as  this 
approach  considerably  simplifies  the  analysis.  In  Part  IX  a  solenoidal  vector 
potential,  for  which  the  divergence  is  zero,  is  derived  by  separation  of 
variables.  Formulas  for  the  vector  displacement  and  its  derivative  remain 
unchanged,  but  a  simpler  formula  for  the  rotation  vector  is  obtained.  The 
usual  approach  by  means  of  the  Newtonian  potential  is  not  appropriate  since 
the  hollow  cylinder  is  a  multiply  connected  region. 


V.  STRAINS,  ROTATION  VECTOR,  AND  STRESSES 


The  strains  are  obtained  from  the  displacements  by  means  of  the  formulas 


e 


r 


9u 

9r 


(68) 


0 


r 


u 


(69) 


e 


z 


9w 

9z 


(70) 
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+ 


3v 

3r 


v 

r 


(71) 


Y 


r0 


3u 

r30 


3u  3w 

Y  “  T"  +  T~ 
rz  3z  3r 


(72) 


3v  3w 
^z0  3z  r30 


(73) 


We  calculate  the  strains  in  terms  of  Bessel  functions  for  Case  I.  The 
remaining  three  cases  can  be  treated  in  a  similar  manner.  Let 


2 

r  e 


2 

r  e 


0 


2 

r  e 

z 


2 


r  Y 


r0 


2 

r  Y 


rz 


2 

r  Y 


0z 


Hi  l  AiV 

6 

H1  \  V2,i 

6 

"i  \ 

6 

H3  \  AiE4,l 
6 

H2  {  AiE5,i 


H4  \  AlE6,i 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


The  j  functions  depend  on  the  radius,  and  can  be  expressed  in  terms  of 
Bessel  functions  of  orders  n  and  n+1. 


“1.1 


r«  2  2  2 
L  s  r 


+  (n  -n)]I  (g  sr) 
n  i 
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6tsrWVr) 


(80) 


(81) 


Rl,2  =  ^Vr2  +  (n2-n)]Kn(81(sr)  +  f^srK^d^sr) 

Rl,3  =  (n2-n)In(P2sr)  +  02srIn+1(P2sr)  (82) 

R1>A  =  (n2-a)Kn(02sr)  -  02srKn+1(02sr)  (83) 

Rl,5  =  “  82s2r2][n(32sr)  +  (n+1>srIn+1(e2sr)  (84) 

Rl,6  =  P2s2r2Kn(e2Sr)  +  (n+1)srKn+1(02sr)  (85) 

R2,l  =  "  (n2“n)In(S1sr)  +  pjsrl^^pjsr)  (86) 

R2,2  =  “  (n2~n)Kn(81sr)  -  ^srK^^sr)  (87) 

2 

R2,3  =  “  (n  -n)In(02sr^  “  P2nSrIn+1^2Sr^  (88) 

R2,4  =  "  (n2_n)Kn(02sr)  +  82nsrKn+1(P2sr)  (89) 

R2,5  =  “  (n+1)srln+1(82sr)  (90) 

R2,6  =  "  (n+1)srKn+1(82sr)  (91) 

r3>1  ■  "  S2r2In(8lSr)  (92) 
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(93) 


R-,  9  =  -  s2r2K  (g.sr) 

J  yJL  n  I 


R3,3  =  0 


(94) 


R3,4  =  0 


(95) 


R_  =  g  s2r2I  (g~sr) 

J  ,3  l  n  l 


(96) 


R9  r  -  ~  39s2r2K  (g.nsr) 

j  ,b  Z  n  1 


(97) 


R^  j  =  -  2(n  -n)In(g1sr)  -  gj,nsrln+1  (g^r) 


(98) 


R,  9  =  -  2(n  -n)K  (g.sr)  +  g  nsrK  (g  sr) 

4  ,2  n  1  1  n+1  I 


(99) 


R^  3  =  -  [2(n2-n)  +  g2s2r2 ]In(g2sr)  +  2g2srln+1 (g2sr) 


(100) 


R^  ^  =  -  [2(n2-n)  +  g2s2r2 ]Kn(g2sr)  -  2g2srln+1 (g2sr) 


(101) 


R^  5  =  -  g2s2r2In(g2sr)  +  2(n  +  l)srln+1 (g2sr) 


(102) 


R4  6  =  ^2s^r^Kn^2Sr^  +  2^n  +  1)srKn+i(^2sr^ 


(103) 


2  2 

R-  ,  =  nsrl  (g.sr)  +  g.s  r  I  (g  sr) 

5,1  n  1  1  n+l  1 


(104) 
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R5j2  =  nsrKn(0lSr)  -  PjsVl^CBjSr)  (105) 

R5>3  =  nsrln(32sr)  (106) 

R5j4  =  asrKn(32sr)  (107) 

R5  5  =  -  32nsrln(32sr)  +  (322  +  l)s2r2In+1(32sr)  (108) 

R5>6  =  32nsrKn(32sr)  +  (3^  +  l)s2r2Kn+1(32sr)  (109) 

Rt  ,  =  -  2nsrl  (3,sr)  (110) 

D ,  I  n  l 

R,  _  =  -  2ns rK  (3,sr)  (111) 

o,2  n  1 

R6,3  =  ”  nsrln(32sr)  -  32s2r2In+1(32sr)  (112) 

R6>4  =  -  nsrln(82sr)  +  32s2r2Kn+1(32sr )  113) 

R6,5  =  e2nsrln(32sr)  -  s2r2In+1(32sr)  (114) 

R6,6  =  "  e2nsrKn(32sr)  -  s2r2Kn+1(32sr)  (115) 


The  dilatation  e  is  the  sum  of  the  principle  strains 
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e 


(116) 


=  e 

r 


+  e0 


+  e 

z 


We  write 


2 

r  e 


(117) 


where 


2  2  2 
s  r 


V6isr> 


(118) 


a.s^r^K  (3.sr) 
I  n  1 


(119) 


We  note  that 


otj  +  0  as  c  *  0 


Hence 


e  +  0  as  c  ->■  0 


In  general,  therefore,  the  solution  in  terms  of  wave  functions  does  not  yield 
a  valid  solution  in  the  limit  as  the  velocity  of  travel  approaches  zero.  Pure 
torsion  is  an  exception,  as  only  shear  strains  are  involved  and  the  dilatation 
is  zero  under  both  static  and  dynamic  loading  conditions. 
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The  rotations  are  given  by  the  formulas 


2w 

z 


l_  _8u.  9v_ 

r  80  8r 


(120) 


0  8u  8w 

8z  3r 


(121) 


2w 


1  3rv 

r  8r 


J_  3u_ 
r  80 


(122) 


We  set 


0  2 
2r  w 


10 

H3  jj  AiR8,i 


(123) 


„  2 
2r  w„ 


10 

H2  l  \l,  t 

z  g  1  *»x 


(124) 


2r2w 


10 


H/  l  A-R,n  • 
4  g  i  10,i 


(125) 


for  Case  I. 


E8,3  "  nstIn(62sr>  +  B2s2r2ln+l(62sr) 


(126) 


R8,4  ‘  nstKn(82sr)  '  82s2r2Kn+l(l32sr) 


(127) 


E8,5  ■  «2"srln(82sr)  +  S2r2In+1(82sr) 


(128) 
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(129) 


R8,6  =  "  32nsrKn(32sr)  +  s2r2Kn+l(B2sr) 

Rg  3  =  nsrln(82sr)  (130) 

Rg>4  =  nsrKn(82sr)  (131) 

R9,5  =  32nsrln(82sr)  -  c^I^^sr)  (132) 

R9  6  =  "  32nsrKn^2sr^  ~  a2s2Kn+1^2sr^  (133) 

E10,3  "  '  S282r2In(62sr)  (134) 

E10,4  --  VZr2,Cn(62sr)  (135) 

E10,5  "  'i”1."!")  <136) 

E10,6  ■  -  <137) 

The  stresses  are  derived  from  the  strains  by  means  of  generalized  Hooke's 

law. 

or  =  Xe  +  2per  (138) 

o0  =  Xe  +  2ye0  (139) 
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a 

z 

Xe  +  2pe 

z 

(140) 

Tr0 

■  ^re 

(141) 

T 

rz 

■  v\z 

(142) 

Y0z 

*  ”6. 

(143) 

For  case  I 

we  write 

2 

r  o 

r 

6 

■  H‘  1 

AiSl,i 

(144) 

2 

r  °e 

6 

■  5 

AiS2,i 

(145) 

2 

r  a 

z 

6 

’  H‘  i 

AiS3,i 

(146) 

2 

r  Tr0 

6 

=  H3  X 

1 

AiS4,i 

(147) 

2 

r  x 

rz 

6 

'  "*  1 

AiS5,i 

(148) 

2 

r  T0z 

6 

=  h  £ 
4  * 

AiS6,i 

(149) 

Then 
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S  =  XR  +  2yR  ,  j  =  1,2,3;  i  =  1,3,5 

J  >  X  '  ,  A  J)1- 

S.  .  =  XR  +  2yR.  .,  j  =  1,2,3;  i  -  2,4,6 

J*1  '  >2  J>i 

S.  .  =  yR.  .,  j  =  4,5,6;  i  =  1,2, 3, 4, 5, 6 

J  >  1  J  >  1 


(150a) 


(150b) 


(151) 


Stresses  and  strains  for  Cases  II,  III,  and  IV  can  be  obtained  in  a  similar 

manner.  The  R.  .  functions  must  be  calculated  for  each  case  to  determine  the 
1  >  J 

correct  algebraic  signs  in  Eqs.  (80)  -  (115)  and  Eqs .  (118)  -  (119).  A  similar 
remark  applies  to  the  formulas  for  the  rotations. 

VI.  RESPONSE  TO  AN  INTERNAL  TRAVELING  PRESSURE  PULSE 

First,  we  consider  a  sinusoidal  pressure  pulse  traveling  with  velocity  c. 

The  boundary  conditions  are 


o ^  =  a cos(n9)sin[s(z-ct) ] ,  r  «  a 


T 


r0 


0,  r  =  a 


x 

rz 


0,  r  =  a 


x  =  0,  r  =  b 

r 


T 

re 


0,  r  =  b 


T 

rz 


0,  r  =  b 


(152) 


(153) 


(154) 


(155) 


(156) 


(157) 
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These  boundary  conditions  lead  to  a  set  of  six  linear  equations  for  the 
coefficients • 


6 

l 

i 

6 

l 

1 

6 

l 

1 

6 

l 

1 

6 

l 

1 

6 

I 

1 


A.  S  . 

l  l,i 


A.  S.  . 
l  4,i 


A.  S  . 
l  5,1 


A.  S.  , 

l  1,1 


A.  S.  , 
l  4,i 


A.  S  . 
l  5,i 


(a) 

(a) 

(a) 

(b) 

(b) 

(b) 


2 

a  a 

o 


0 


0 


0 


0 


0 


058) 


(159) 


(160) 


(161) 


(162) 


063) 


These  equations  can  be  solved  for  the  A^  provided  c  is  greater  than  zero  and 
less  than  the  velocity  of  the  shear  wave,  according  to  Eq.  (20).  It  is 
convenient  to  write  the  final  results  in  terms  of  determinants,  using  Cramer's 
rule  and  then  combining  the  separate  terms.  Let  D(s)  be  the  determinant  for 
Eqs .  (158)  -  (163): 
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1 , l(a) 

sl>2(a) 

Sl,3(a) 

sl>4(a) 

sI>5(a) 

Sl,6<a> 

;4,l(a) 

S4,2(a) 

s4j3(a) 

S4;4(a) 

S4_5(a) 

S4,6(a) 

5,  l(a) 

S5,2(a) 

S5,3(a) 

S5>4(a) 

S5>5(a) 

%,6(a) 

Sl,2(b) 

slj3(h) 

Sl,4<b> 

Sl,5<b> 

sl,6<b> 

4,l(b) 

s4>2(b) 

s4j3(h) 

S4>4(b) 

s4>5(b) 

s4,6<b> 

5,l<b> 

s5j2(h) 

S5 , 3(b) 

s5>4(b) 

s5>5(b) 

S5,6<b> 

(164) 


The  remaining  determinates  use  the  appropriate  R.  •  or  S.  .  for  the  first  row; 
the  remaining  five  rows  are  identical  with  the  corresponding  rows  of  D(s).  We 
have, for  instance, 


Sr  = 


^  i ( r )  2^r ^  ^(r)  Si^^(r)  Sj^tj(r)  $i  ^(r) 


1,6s 


r0 


rz 


S4,l(r)  S4,2(r)  S4,3(r)  S4,4(r)  S4,5(r)  S4,6(r) 


Rl,l^r^  R1 . R1 , 3^r  ^  R1 ,4^r^  R1 , 5^r ^  R1 , 6^r ^ 


(165) 


(166) 


(167) 


Determinates  for  the  remaining  stresses  and  strains  can  be  written  down  by 
inspection.  The  stresses  for  Case  I  are  given  by 


2 

r  a 

r 


H1Sr/D(s) 


(168) 


2 

r  a 


0 


H1SQ/D(s) 


069) 
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H, S  /D(s) 
1  z 


(170) 


2 

r  a 

z 


2 

r  t 


r9 


2 

r  t 


rz 


2 


r  x 


0z 


H3Sr0/D(s) 

H  S  /D(s) 
z  rz 


(171) 


(172) 


(173) 


Formulas  for  the  strains  are  obtained  in  a  similar  manner. 


2 

r  e 

r 


2 

r  £e 


2 

r 


t. 

z 


2 

r  Y 


rz 


2 

r  Y 


6z 


HjRVdCs) 


HjRg/DCs) 


H.R  /D(s) 
I  z 


H3Rr0/D(s) 


H2Rrz/D(s) 


H4R8z/D(s) 


(174) 


(175) 


(176) 


(177) 


(178) 


(179) 


We  can  also  solve  problems  specified  by  the  boundary  conditions 


29 


r 


a 


(180) 


=  Tocos(n0)sin[s(z-ct) ] ,  r  =  a 


(181) 


t 

rz 


=  0  , 


r  =  a 


(182) 


and 


a  =  0  ,  r  -  a 

r 


(183) 


T 


r0 


0  , 


r  =  a 


(184) 


t  =  t  cos(n0)sin[s(z-ct )z] 
rz  o 


(185) 


It  is  assumed  the  outside  surface  is  stress-free  according  to  Eqs ,  (155)- 
(157)*  The  and  functions  must  be  derived  for  each  case. 

Formulas  for  the  response  to  a  traveling  step  function  of  pressure  are 
obtained  by  writing  the  stresses  and  strains  in  terms  of  Fourier  integrals.  In 
place  of  Eq.  (152)  we  write 


a 

r 


a  cos(n0) 
o 


[Vi 


J 

O 


sin [s(z-ct ) ] 
sD(s) 


ds  ] 


r  =  a 


(186) 


Then 


°r  *  °o  cos(n^)  »  z>ct 


(187) 
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°v  "  V2  aQ  cos(n0)  ,  z  =  ct  (188) 


-  0  ,  z<ct  (189) 


The  terra 


a  =  Vo  a  cos(n0)  ,  r  =  a  (190) 

r  ^  o 


leads  to  elementary  solutions  obtained  by  evaluating  the  corresponding  Fourier 
integral  at  the  origin.  The  number  of  terras  depends  in  part  on  the  multiplicity 
of  the  zero  at  the  origin  of  the  complex  s  plane. * 


2 

r  a 

r 


0  00 

—  cos(n0)  / 

TT  J 

o 


sin[s(z-ct)]  j 
sD(s) 


(191) 


leads  to  Fourier  integrals  involving  functions  occurring  on  the  right-hand  side 
of  Eqs  •  (173)  -  (184),  We  have 


0 

r 


—  cos(n0) 

TT 


/ 

O 


sin[s(z-ct) ]  „ 

- — r -  b  dS 

sD(s)  r 


(192) 


0 


0 


1 

TT 


cos(n0) 


1  s» ds 

;  sD(s)  0 


(193) 


0 


z 


1 

TT 


cos (n0 ) 


/ 

o 


HslstegXI.  s  ds 
sD(s)  z 


(194) 


A Appendix  7T.  These  elementary  solutions  can  be  expressed  in  terms  of  zonal 
harmonies . 
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00 

1  j  /  Ck  \  f  sin[s(z-ct) ]  0  , 

r0  tt  v  J  sD(s)  r9 


(195) 


rz  7T 


cos(n6)  / 
o 


-■n^v-T)1  s  ds 

sD(s)  rz 


(196) 


u  0z 


00 

1  •  /  r  cos[s(z-ct) ]  ^ 

-  —  Sin(n0)  J  - - — r  SQ  ds 

it  J  sD(s )  0z 


(197) 


The  corresponding  formulas  for  the  strains  are 


r  it 


1  /  a\  r  sin[s(z-ct)  ]  _  - 

Z  cos(n0)  J  - -  r  /J  R  ds 

w  J  sD(s )  r 


(198) 


oo 

/  a \  r  sin[s(z-ct)]  , 

=  -  cos(n9)  /  - sD(s)  R0  ds 


(199) 


00 

1  ,  QN  r  sin[s(z-ct)]  _ 

=  —  cos  (n0 )  J  - ■  N  R  ds 

TT  J  sD(s)  Z 


(200) 


r0 


1  .  /  Q\  r  sm[s(z-c)]  _  , 

—  sin  n0  I  - ■  7-7  v  '  R  Q  ds 

it  J  sD(s)  r0 

o 


(201) 


rz 


2 

—  cos 

TT 


j  x  ds 

J  sD(s)  rz 


(202) 


2 

t  =  —  sin 

0Z  TT 


/  —  o \  f  cos[s(z-ct)] 

(n0)  J  - Sd (T) — R0z  ds 

o 


(203) 


The  total  value  of  the  stresses  and  strains  is  obtained  by  adding  solutions  to 
the  problem  defined  by  Eq.  (174). 

Solutions  to  problems  involving  other  boundary  conditions  may  be  obtained 
in  a  similar  manner.  In  certain  cases  the  Fourier  integrals  will  be 
divergent.  In  these  cases  the  limits  of  integration  should  be  taken 
between  -00  and  +  00 ,  and  the  Cauchy  principal  value  calculated.  The  factor 
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2/  7T  must  be  replaced  by  1/  tt  in  these  cases. 


VII.  HARMONIC  AND  BIHARMONIC  FUNCTIONS  AS  LIMITS  OF  WAVE  FUNCTIONS 

As  the  velocity  of  travel  approaches  zero,  the  wave  functions  tend  to 
harmonic  functions  in  the  limit,  as  shown  below.  However,  in  general,  a  static 
problem  requires  biharraonic  functions  as  well  as  harmonic  functions  for  a 
complete  solution**  A  second  limiting  process  involving  an  indeterminate 
quotient  is  required  to  extract  the  biharraonic  functions  from  the  wave 
functions . 

We  recall  the  basic  equations  of  Part  II. 


u  =  V<j)  4-  V  x  ij; 


2  ii 

C1  V  ♦  =  2 
3t 


2  _2  ?  3  > 

c9  v  ♦  =  — £ 

at 


(5) 


(6) 


(7) 


We  find  from  Eqs .  (19),  (29),  (30),  (31),  (32),  and  (23)  that 


3\J; 


3t 


2 


2  2a 

c  s  9 


2  2* 
C  S  t 


(204) 


(205) 


On  letting  c  +  0  we  find 


WvoFUems  involving  pure  shear,  such  as  torsion  of  an  axisymmetvic  solid ,  are 
exceptions.  A  scalar  biharmonic  function  is  required  in  problems  involving  a 
change  in  volume,  since  dilatation  derived  from  a  scalar  harmonic  function  is 
zero . 
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0 


(206) 


V2if ,  =  0 


where 


<|>  =  [A. I  (sr)  +  A.K  (sr)]cos(n0)sin(sz) 

In  l  n 

A  <<\ 

ip  =  e  [ACI  . ,  (sr)  +  A,K  (sr)]cos(a0  )cos(sz) 
r  5  n+1  o  n+l 


-  eQ  A^I  ^  (sr )  +  (sr )  cos(n0)  sin(sz) 


+  e  [A0 I  (sr)  +  Ay K  (sr ) ] sin(n0 )sin(sz) 
z  3  n  4  n 


The  corresponding  displacements  are  found  by  setting  c  ■  0,  3^  =  1,  3 
Eqs.  (64)  -  (66)*  We  can  eliminate  two  of  the  A^  coefficients  in  the 
writing 


=  Ax  + 


A8  =  A2  +  A6 


A 


9 


A 


5 


A6 


(207) 


(208) 

(209) 


=  1  in 
result  by 


(210) 

(211) 

(212) 


(213) 


Then 


ru  =  [A^nl  (sr)  +  A0nk  (sr)  +  A^rl  ,.(sr)  -  A,_srK  , ,  (sr)]cos(n9)cos(sz) 

/  n  on  y  n+i  J.U  n+l 

(214) 


rv  =  [-A7nln(sr)  -  AgnKn(sr)  -  Agsrln+1(sr)  +  A10srKn+1(sr) ]sin(n9)sin(sz) 

(215) 


rw  =  [A?  -  Ag)srln(sr)  -  (Ag  +  A^g)srKn(sr) ]cos(n9)cos(sz)  (216) 


We  note  that  the  six  independent  solutions  for  the  displacements  for  the  wave 
equation  reduce  to  four  when  the  velocity  of  travel  becomes  zero.  Moreover,  the 
displacements  vanish  if 


A±  =  0,  i  =  7,8,9,10  (217) 


To  obtain  additional  solutions,  let 


A.  =  1/(61  -  B2)  ,  i  =  1,2  (218) 


in  Eqs .  (217),  (218),  and  (219).  We  define  the  functions 


Fl  =  (nIn(6lSr) 

F  =  [nK  (8  sr) 
l  n  1 


nln(82sr)]/(81  -  $2) 
nKn(82sr)]/(81  -  82) 


(219) 


(220) 
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F3  =  [eisrln+1(sr)  -  srln+1(sr)(e2  -  sr)]/(B1  -  32>  (221) 

=  [gjSrK^jCe^r)  -  srK^^sr)]/^  -  32>  (222) 

F5  =  [62srln+1(02sr)  -  srln+1(32sr) ]/(01  -  32)  (223) 

F6  ‘  162srKn+l(S28r)  -  srK„+l(62sr)1/(6l  -  S2>  (224) 

F  =  [srl  (3,sr)  -  39srl  (39sr)]/(3,  -  39)  (225) 

/  nl  l  n  2  1  l 

F  =  [srK  (3,sr)  -  3~srK  (30sr)]/(3,  -  30)  (226) 

o  nl  l  Vi  l  12 


On  referring  to  Eqs.  (64),  (65),  and  (66)  we  find 


ru  «  +  F3^  +  ^2^F2  "  F^) ]cos(n0)sin [s ( z-ct ) ]  (227) 


rv  =  [-  31(F1  +  F  )  +  ^ 2 ( F 2  +  F^)]sin(n0)sin[s( z-ct ) ]  (228) 


rw  =  [3  (Fg  -  62Fg]cos(n0)cos[s(z-ct)]  (229) 


Each  of  the  F^  functions  has  the  indeterminate  form  0/0  when  c  = 

0,0^  =  1,  02  =  lj  this  is  obviously  also  true  for  the  displacements.  The 

limits  of  these  indeterminate  forms  can  be  found  by  the  ordinary  rules  of 
calculus.  When  the  velocity  c  is  small  we  have 
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1 


1 


(230) 


6,  -  V2  a/  , 


-  3, 


V, 


o„ 


and 


lim 

c+0 


1  -  3 
1  -  3 


J, 

2 


h 

X  +  2p 


on  referring  to  Eqs.  (2),  (3),  (22),  (23),  (27),  and  (28).  The  limits  of 
functions  can  be  found  by  using  two  terras  of  Taylor's  series: 


^(Bjsr)  =  In(sr)  +  -  l)srl'(sr)  +  6(pJ  -  l)2 

Kn(BlSr)  =  Kn(sr)  +  (gj  -  l)srK'(sr)  +  5(8j-  l)2 

Vi(ti8r)  ■  wr>  +  <6i  -  i)sri;+1(sr)  +  5<bi  - 1>2 

Kn+l(3lsr)  =  Kn+l(sr)  +  (31  "  Dsrr+1(sr)  +  6(Bj  -  l)2 

Additional  formulas  are  obtained  by  substituting  8^  for  8^« 

As  c  approaches  zero  we  find 


lira  F,  =  nsrl'(sr) 
1  n 


lim  F  =  nsrK' (sr) 
2  n 


(231) 


the  Fj^ 


(232) 


(233) 


(234) 


(235) 


(236) 


(237) 
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(238) 


lim  F, 


t — r —  srl  ,.(sr)  +  srl'(sr) 

A  +  y  n+1  n+1 


lim  F, 


TT7srK„+i(sr)  +  srWsr) 


(239) 


lim  Fr 


A  +  2y  T  T  ,  N 

t— T -  srl  I  (sr) 

A  +  ii  n  n+1 


(240) 


lim  F, 


A  +  2y 


srK  (sr) 


X  +  y  n+1 


(241) 


lim  F_ 


^  srl  (sr)  +  srl'(sr) 
X  +  y  n  n 


(242) 


lim  F 


8 


T  t  j~y  srK (sr)  +  srK'(sr) 
A  +  y  n  n 


(243) 


The  derivatives  may  be  eliminated  by  means  of  Eqs .  (60)  -  (63).  The 
displacements  corresponding  to  these  limits  are 


ru  =  lim  [3j(Fj  +  F^)  +  82(F  — ^ 4) ]cos(n9)s in [ s ( z— e tt ) ]  (244) 

c-H) 


rv  =  lim  [-  8  (F  +  F  )  +  09(F  +f  ) ]sin(n0)sin[s(z-ct)]  (245) 

c-K)  11  J  z  z  0 


rw 


lim  [0  F 
c-K)  ' 


02Fg]cos(n0)cos [s(z-ct)] 


(246) 


The  total  displacements  due  to  static  loading  are  obtained  by  adding  these 
displacements  to  the  displacements  obtained  from  the  harmonic  potentials,  Eqs. 
(217)  -  (219).  We  now  have  a  total  of  six  linearly  independent  solutions  for 
the  displacements.  The  strains  can  be  derived  from  the  displacements  and  the 
stresses  follow  immediately  from  Eqs.  (156)  and  (151).  Thus,  six  independent 
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stress  formulas  are  obtained  for  meeting  six  boundary  conditions  of  the  type 
given  by  Eqs.  (152)  -  (157)  with  c  =  0. 

VIII.  BIHARMONIC  SCALAR  AND  VECTOR  POTENTIALS  FOR  STATIONARY  LOADS 

It  is  possible  to  derive  the  displacements  given  by  Eqs.  (250)  -  (252)  from 
vector  and  scalar  potentials  according  to  Eq.  (5).  Harmonic  scalar  and  vector 
potentials  are  given  by  Eqs.  (211)  and  (212).  We  now  derive  the  additional 
biharmonic  potentials  required  for  Eqs.  (250)  -  (252). 

For  the  scalar  potential,  assume 


X  = 


(247) 


Then 


ix 

3r 


r 


9(j> 

3r 


(247a) 


or 


9x  _  1  324> 

9r  r  392 


(247b) 


since  <{)  is  a  scalar  harmonic  function  it  follows,  after  some  routine  analysis, 
that  ^ 


2 

VZX 


(248) 


y - - - 

A.E.H.  Love 3  A  Treatise  on  the  Mathematical  Theory  of  Elasticity ,  Dover 
Publications 3  New  York 3  1944.  See  pages  274-277 ~  especially  Eqs.  (66)  and  (67). 
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and 


4 

VHX  -  0 


(249) 


Other  scalar  biharmonic  functions  are  given  by  References  8,  9,  and  10. 


2 

x  =  r  ♦  (250) 

X  =  (r  sin0)<f>  (251) 

X  =  (r  cos 0 ) <J>  (252) 

x  ■  (253) 


However,  the  biharmonic  functions  given  by  Eqs.  (250)  -  (251)  seem  to  be  most 
closely  related  to  scalar  potentials  for  wave  motion,  as  they  can  also  be 
derived  from  a  scalar  wave  function  by  a  limiting  process.  We  have,  for  instance, 


8 


1 


(254) 


(i’C’  Fung3  foundations  of  Solid  Mechanics ,  Prentice  Hall,  Inc.,  Enqlewood 
Cliffs,  NJ,  page  208. 

9H.  Neuber,  Theory  of  Notch  Stresses :  Principles  for  Exqct_Stress 
Ccurul^^  David  Taylor  Model  Basin,  Washington,  DC,  Translation  74 
h Jovember  1945.  See  pages  25  and  128. 

10A-S.  Elder,  "Traveling  and  Stationary  Loads  on  the  Half  Space,"  BRL  report 
to  be  published.  See  section  titled  "Biharmonic  Functions  as  Limits  of 
Wave  Functions.  " 


40 


l/2  [AiIn(Blsr)  +  A2Kn(BlSr)] 


(255a) 


4> 


1 


x  [sin  s(z  +  ct)  +  sin  s (z-ct )  ] cos( n0) 


or 


<j>  =  [A  I  (0  sr)  +  A  K  (0  sr)]  sin(sz)cos(sct )cos(n0) 

X  X  ^  x  Z  ^  X 


(255b) 


This  scalar  potential  represents  two  equal  loads  traveling  with  the  same  speed 
in  opposite  directions. 


3<h 

=  V2  [Aj3^rl^(3^sr)  +  A  sj)l  [sin[s(z+ct)3  +  sin[s(z-ct)]cos(n0) 

(256a) 


+  V2  [A1In(31sr)  +  A2Kn(31sr)] 

x  [(z+ct)cos  s(z+ct)  +  (z-ct)cos(sz)-ct) ]cos(n0) 

a<)) 

=  [A^3^rl^  (3^sr)  +  A^^rR/  (3jSr)]  sin(sz)cos(sct)cos(n0) 

+  [A^I^(3jSr)  +  A2K^(3 ^sr ) ] [z  cos(sz)]  cos(sct)cos(n0)  (256b) 


-  [AjI^(3jSr)  +  A2Kn(3^sr) ] [ct  sin  (sz)]  sin(sct)cos(n0) 


On  allowing  c  to  approach  zero,  we  find 


9(J) 

g^—  =  [A^rl^(sr)  +  A2rK^(sr)]sin(sz)cos(n0) 
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(257a) 


+  [A^zl^(sr)  +  A^zK^Csr) ]zcos(sz)cos(n0) 


We  also  have 


X  = 


-5-  [A.  I 
3r  L  1  n 


(sr)  + 


A^K^Csir)  ]  sin(sz)cos(n0  ) 


or 


X  =  [A^srl^(sr)  +  A^srK^Csr ) ]sin(sz)cos(n0) 


which  yields  two  terms  in  Eq .  (256b)  except  for  a  factor  s.  The 
terms  of  (256b)  correspond  to  Eq.  (253). 

The  biharmonic  vector  potential  can  be  derived  in  a  similar 
assume 


3r 


where 


0)  = 


e  w 
r  r 


+ 


e0U)0 


e  a) 
z  z 


Then 


=  r 


3r 


A  2 


(257b) 


remaining  two 

manner.  We 


(258) 


(259) 


(260) 


r 


(261) 


0)„ 


If  *8 


0) 


r  sF^z 


(262) 


We  find  as  before 


a) 

dr  r 


=  r 


3r 


2  *r  +  3?  *r 


(263) 


with  similar  equations  for  0)^  and  The  second  derivatives  are  eliminated  by 
using  the  three  equations  which  result  when  the  derivatives  with  respect  to  time 
are  set  equal  to  zero  in  Eqs.  (11),  (12),  and  (13).  On  referring  to  Eq.  (8)  we 
find 


—  U) 

dr  r 


3r 


1  .  ,23 

—  ^  +  —  -r— 

r  r  r  dr  r 


I  1 

r 


-  r 


7*6 


2  3 

r  3r 


1  3 

r 


2  yr 


3z 

„2 


2 


30 


2  *0 


-  r 


3z 


2  *0 


(264) 


(265) 


and 


-5—  0) 
d  2  Z 


1 

r  30 


2  *z 


-  r 


3z 


2  +z 


(266) 


We  finally  obtain 


„2  1  1  3 

7  “t  ‘  T“r  +  “2  36  “e 

r  r 


3z 


2  r 


2 


(267) 


1 


(268) 


v2  *c 


2  “e 


+ 


ae  “e 


e 


and 


V 


2 


0) 

z 


2 


(269) 


Consequently, 


V 


2 


0) 


(270) 


and 


V 


4 


(271) 


The  biharmonic  scalar  and  vector  potentials  are  not  independent,  but  are 
coupled  through  Navier's  equation. 


(X  +  2p)V<J)  4-  yVxoj  =  0 


(272) 


or 


[(X  +  2y)  V<j>  +  uVxiJ]  =  0  (273) 

3z 


This  equation  is  satisfied  if 


(X  +  2\i )  Vtf)  +  nVxij> 


0 


(274) 


or,  in  scalar  form, 

(X  +  2u)  ft  +  „  r  fj  *s  -  fj  *e  -  0  (275) 

<X  +  2K)  r  ff  +  „  fj  *r  -  fj-  -  °  (276) 

(X  +  2»)  ff  +  »  f  ff  *e  -  j  fj  *r  -  0  (277) 

We  use  Case  I  as  an  example • 


♦  =  [A  I  (sr)  +  A  K  (sr)]cos(n0)sin(sz)  (278) 
in  in 

“  [A5In+i(sr)  +  A6Kn+1(sr)]sin(n0)cos(sz)  (279) 
3  “  [A5In+i(sr)  +  A6Kn+1(sr)]cos(n0)cos(sz)  (280) 
=  [A^I^(sr)  +  A^K^(sr) ]sin(n0)sin(sz)  (281) 


Relations  among  the  A^  components  are  obtained  by  substituting  these  scalar 
components  into  Eqs.  (11)  -  (14)  and  carrying  out  the  indicated  analysis. 


We  finally  obtain 


X 


(282) 


[3nsrr(sr)  +  812(sr)]cos(n9)sin(sz) 

=  “  [P^srl^sr)  -  &12srK'+1(sr)]sin(n0)cos(sz)  (283) 

^0  =  “  t^nsrIn+i^sr^  +  ^12SrKn+l^Sr)  ]cos(n9)cos(sz)  (284) 

*z  =  “■  [-0usrr(sr)  +  812srr(sr)]sin(n0)sin(sz)  (285) 

where  3^  and  3^  are  related  linearly  to  the  coefficients.  The  derivatives 
may  be  eliminated  by  using  standard  formulas,  so  that  the  final  result  can  be 
expressed  in  terms  of  trigonometric  factors  and  Bessel  functions  of  order  n  and 
n+L  • 


Thus,  we  have  a  complete  solution  for  static  loading  in  terras  of  scalar  and 
vector  potentials.  By  contrast,  Love's  strain  function  uses  a  single  biharraonic 
scalar  potential  based  on  the  Galerkin  vector.  The  displacements  in  our 
analysis  are  given  in  terms  of  first  partial  derivatives  of  the  potential 
functions,  consistent  with  the  Green-Lame  formulation  of  wave  motion,  whereas 
love's  strain  function  expresses  displacements  as  second  derivatives  of  a  scalar 
potential.  The  present  analysis  illustrates  a  logical  connection  between  static 
and  dynamic  problems  for  the  elastic  deformation  of  a  hollow  cylinder. 


IX.  HELMHOLTZ  THEORY 


Helmholtz's  theorem  concerning  resolution  of  a  vector  consists  of  two  parts: 


u  =  V<J>  4-  Vx^ 


(286) 


V  •  X  - 


0 


(287) 


The  vector  potential  of  Eq.  (5),  as  derived  in  the  previous  analysis,  is 
generally  not  solenoidal,  as  required  by  Eq,  (287).  A  solenoidal  vector  is 
obtained  by  a  gauge  transformation,  as  outlined  by  McQuistan.^  In  the  course 
of  McQuistan's  analysis  a  Newtonian  potential  is  used  to  solve  Poisson's  partial 
differential  equation,  and  the  proof  is  restricted  to  simply  connected 
regions.  Unfortunately,  a  hollow  cylinder  is  multiply  connected,  and  this  proof 
does  not  apply.  We  solve  Poisson's  equation  by  separation  of  variables.  We 
have  three  cases  to  consider: 


1)  i|)  is  a  harmonic  vector  potential 


A 

2)  iji  is  a  biharmonic  vector  potential 


3)  iJj  is  a  vector  wave  function. 


In  each  case  we  assume 


x  =  i})  +  Vn 


(288) 


where  n  is  a  scalar.  Then 


V  • 


X 


v  •  v  n 


(289) 


which,  in  view  of  Eq.  (287),  reduces  to 


U 


R.B.  MeQuistan,  Scalar.  and  Vector.  Fields,  A  Vhxjsical  Interpretation^, 
Wiley  and  Sons ,  Inc. ,  New  York,  1965.  See  pages  256-264. 


John 
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(290) 


V 


2 


n 


v  • 


* 


Since  the  vector  potential  i|>  has  been  calculated  previously,  Eq.  (290)  is 
Poisson's  equation  for  the  scalar  t|  .  In  each  case  we  have 


<J>  =  -  (ri|>  )  +  - 

y  r  3r  r  r 


W+9  + 


where 


e  +  e  _  +  e  ip 

r  r  0  0  z  z 


A  A  A 

and  e^,  e^,  and  e z  are  the  unit  vectors  in  cylindrical  coordinates.  For 
we  have 


=  [A5In+1(sr)  +  A6Kn+1(sr)]sin(n0)cos(sz) 


=  -  [A5In+1(sr)  +  A6Kn+1(sr)]cos(n0)sin(sz) 


ib  =  [ Ac I  (sr)  +  A.K  (sr)  ]sin(n0)sin(sz) 

Yz  5  n  4  n 


After  some  routine  analysis  we  find 


V  •  =  [Ag  +  A^)sln(sr)  +  (A^-A^)sKn(sr)]sin(n0)cos(sz) 


(291) 


(292) 


Case  1 


(293) 


(294) 


(295) 


(296) 
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We  use  the  semi-inverse  method  to  solve  Eq.  (290).  This  was  a  lucky  guess  which 
saved  considerable  analysis .  Assume 


n  =  C  [A_  +  A  )I'(sr)  +  (A.  -  A  )K'(sr)]sin(n0)cos(sz)  (297) 

J  j  n  4  on 

where  C  is  an  unknown  constant • 

Then 

V2n  =  2  C  s2[(A0  +  AC)I  (sr)  +  (A,  -  A.)K  (sr)  ]  sin(n9)cos(sz)  (298) 

Jon  4  o  n 

On  referring  to  Eqs •  (290)  and  (291)  we  see  that 

c  =  -  V2S  (299) 

and,  consequently, 

n  *  -  Vo  [(A0  +  Ac)I'(sr)  +  (A,  -  A,)K  (sr)  ]  sin(n0)cos(sz)  (300) 

^  J  j  n  4  o  n 

The  biharmonic  vector,  Case  2,  may  be  treated  in  a  similar  manner.  We  have 

^r  =  [01srr+1(sr)  -  32srK^+1(sr)]sin(n0)cos(sz)  (301) 

=  ~  [0jSrr+j(sr)  -  32srK^+1(sr)]sin(n0)cos(sz)  (302) 


49 


[-  3,sr'(sr)  +  3„srK' (sr) ]sin(n0)sin(sz) 

In  L  n 


(303) 


♦ 


z 


We  f ind 


V*i|>  =  [3^sln(sr)  +  32sKn(sr)]sin(n0)cos(sz) 


(304) 


Hence, 


v2n  = 


[3^sln(sr)  +  32SKn(sr)]sin(n0)cos(sz) 


(305) 


On  proceding  as  in  the  previous  case,  we  find 


n 


2  I6lrln<sr> 


+  ^2rKn(sr ) ]sin(n^)cos(sz) 


(306) 


This  result  is  well  known.  If  the  right-hand  side  of  Poisson's  equation  is  a 
harmonic  function,  biharraonic  functions  can  be  used  to  solve  Poisson's 
Equation.  For  instance,  if 


v2n  - 


<f>> 


(307) 


and 


n 


(308) 


then 
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(309) 


n2  9  <J> 

V  n  =  - ir 

9z 


and  we  can  readily  prove  that  n  is  biharmonic, 


For  the  wave  function,  Case  3,  we  have 


=  [ A5 In_(_2 ( ^ 2Sr )  +  A6Kn+1(&2sr^Sin^n0^COS  s(z-ct^  (310) 


<^Q  =  -  [A5In+1(32sr)  +  A6Kn(62sr)]cos(n9)cos[s(z-ct)] 


(311) 


ij>  =  [A  I  (3„sr)  +  A  K  (3  sr)]sin(n9)sin[s(z-ct)] 

z  j  n  Z  4  n  z 


(312) 


where  \J>r,  and  are  components  of  the  vector  potential.  We  have  shown  that 


V  •  -  [(A^  +  32A5)sIn(32sr)  +  (A^  -  32A6)sKn( 32sr )  ]  sin(n9)cos  [s(z-ct )  ] 

(67) 


Assume 


n  =  C  [(A3  +  32A5)sIn(32sr)  +  (A^  -  32A6)sKn(32sr)]sin(n9)cos[s(z-ct)] 


(313) 


where  the  constant  C  is  unknown.  Then 
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V2n  =  -  {C  (e2  -  l)s2  (A3  +  32A5)sIn(S2sr)  +  (A4  -  e2A6)Kn(B2sr) ] }  (314) 


x  {sin(n9 )cos [s(z-ct)z] } 


But 


c^2  =  1  -  322  so  that  (315) 

9  o 

C  =  -  l/s  a2 


and 


n 


1  ■-  [  (A_  +  30AC)I  (30sr)  +  (A.  -  30A,)K  (3~sr ) ] sin(n9 )cos ] [ s(z-ct) ] 
2  3  25n2  4  2  o  n  2 

sct2  (316) 


This  result  is  not  valid  for  =  0,  corresponding  to  zero  velocity  of  travel. 

A 

Next  we  recompute  the  vector  displacement  by  substituting  x  from  Eq.  (288) 
into  Eq.  (286)*  We  find 

A  A 

U  -  V4>  +  Vxij,  +  V.Vn  (317) 


The  third  term  on  the  right-hand  side  of  this  equation  is  zero,  so  we  recover 
Eq.  (5)  in  the  formula  for  the  vector  displacement. 


u  -  V*  +  Vx\Jj  (5) 


Formulas  for  the  stresses  and  strain,  which  involve  partial  derivatives 
of  u,  remain  unchanged. 


The  vector  potential  is  given  by 


0) 


=  V2  Vxu 


(318) 


or 


0) 


=  V2  V  x  V<j>  +  V2  V  x  V  x  i|/ 


(319) 


But 


V  x  =  0 


(320) 


and 


V  x  V  x  ip  =  V  (V»\|>)  -  V  x 


(321) 


Hence. 


1  2" 
~  \  v  * 


(322) 


An  analogous  formula  is  found  for  the  dilatation 


e  =  V»u 


(323) 
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e 


V  •  V<j)  +  V*Vx^ 


(324) 


But 


V  •  Vx^i  =  0 

So  that 


e 


(325) 


(326) 


The  results  obtained  by  using  x  in  place  of  \|>  conform  to  all  requirements 
of  Helmholtz  theory  provided  the  z  axis  is  excluded  from  the  region  under 

consideration  and  the  Bessel  functions  are  defined  to  be  single  valued.  The 

scalar  potential  n  is  determined  so  that  Eq.  (287)  is  satisfied  and  x  is 

A 

solenoidal.  The  potential  x  is  a  single  valued  function  of  0  since 

sin(n  0)and  cos(n  0)  are  periodic,  and  are  obviously  single  valued  functions 

of  z  and  t  if  s  is  real  and  positive. 

If  s  is  complex,  the  complex  s  plane  must  be  cut  along  the  negative  real 
axis  to  make  the  modified  Bessel  functions  of  the  second  kind  single  valued. 
This  cut  makes  the  logarithms  occurring  in  the  definition  of  these  functions 
single  valued.  In  analysis  of  equilibrium  problems,  the  final  equilibrium 
solution  does  not  involve  logarithms  of  s  and  becomes  single  valued;  this  is 
also  true  of  dynamic  problems  involving  Bessel  functions  of  the  second  kind. 
The  proof  requires  additional  analysis  which  is  given  in  Appendix  A. 

X.  DISCUSSION  AND  CONCLUSIONS 

In  this  paper  we  have  formulated  the  equations  governing  elastic  strains 
in  a  hollow  cylinder  due  to  stationary  loads  and  loads  moving  with  constant 
subsonic  velocity.  In  addition,  scalar  and  vector  potentials  are  derived  for 
stationary  loads,  showing  the  connection  between  static  loading  and  the  Green- 
Lame  formulation  of  wave  motion.  Six  linearly  independent  solutions  are 
obtained  for  the  moving  solutions,  corresponding  to  the  six  boundary 
conditions  on  the  cylindrical  surfaces.  However,  when  the  velocity  of  travel 
is  set  equal  to  zero,  two  solutions  are  lost  by  confluence  of  the  solutions, 
and  only  four  linearly  independent  solutions  remain.  Two  additional  solutions 
are  obtained  by  a  limiting  method,  so  a  total  of  six  linearly  independent 
solutions  is  available.  Biharmonic  scalar  and  vector  potentials  are  also 
derived  from  which  these  two  additional  solutions  can  be  calculated. 

The  distribution  of  characteristic  roots  is  not  yet  determined.  We 
speculate  the  roots  in  the  first  quadrant  of  the  complex  s  plane  lie  between 
the  imaginary  axis*  and  a  smooth  curve  passing  through  the  complex  eigenvalues 


*Free  vibrations  of  a  hollow  cylinder  lead  to  real  values  of  u)  the  circular 
frequency,  and  real  eigenvalues  when  the  eigenvalues  are  expressed  in  Bessel 
functions  of  the  first  kind,  provided  a)  is  not  too  large.  Pure  imaginary 
eigenvalues  arise  when  modified  Bessel  functions  are  used. 


for  static  loading.  To  verify  this,  the  characteristic  roots  for  a  traveling 
load  on  a  solid  rod  will  be  calculated  first,  as  the  calculations  should  be 
relatively  simple,  before  investigating  the  eigenvalues  for  a  hollow  cylinder. 

The  planned  programming  will  follow  the  general  pattern  developed  for 
static  loading  of  a  hollow  cylinder,  but  simplified  and  streamlined  to 
expedite  the  calculations.  Each  type  of  loading  will  be  considered 
separately . 

Torsional  loading  will  be  considered  first  due  to  the  simplicity  of  the 
analysis.  Axisymmetric  loading  will  be  considered  next,  as  an  extension  of 
our  thick -walled  cylinder  analysis. 

The  analysis  in  Parts  I  -VIII  is  based  on  the  work  of  Herrmann  and  his 
associates,  in  which  the  vector  potential  is  not  required  to  be 
solenoidal.  In  order  to  conform  to  the  classical  Helmholtz  theory,  a  new 

A 

vector  solenoidal  potential  x  is  derived.  The  previously  derived  formulas  for 
the  displacements,  strains,  and  stresses  are  not  changed.  The  formula  for  the 
vector  rotation  is  simplified  and  becomes  the  vector  counterpart  of  the 
scalar  formula  for  the  dilatation. 
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APPENDIX  A 


QUASI-STATIC  SOLUTIONS  OF  THE  WAVE  EQUATIONS 


Stresses  given  in  Part  V  of  this  report  have  the  generic  form 


a 


G( 0 , ‘ r ,  z)  +  F(0,  r,  z) 


(Al) 


where  G(  0,  r,  z)  is  the  limiting  form  of  the  stress  at  a  great  distance  from 
the  discontinuity  in  loading  and  F(  0,  r,  z)  is  a  Fourier  integral  giving  the 
local  effects  of  this  discontinuity.  Both  solutions  can  be  obtained  from 
scalar  and  vector  wave  functions.  However,  G(  0,  r,  z)  is  of  a  simpler  type, 
and  can  be  obtained  from  a  scalar  harmonic  function  in  the  independent 
variables  0,  r,  (z-ct)/0^  and  a  vector  harmonic  function  in  the 
variables  0,  r,  (z-ct)/02-  These  solutions  can  be  expanded  in  power  series 
in  the  variable  (z-ct);  logarithmic  solutions  may  also  occur.  Exponentially 
decaying  terms  do  not  occur,  so  the  stresses  given  by  G(  0,  r,  z)  persist  at 
considerable  distances. 

The  scalar  wave  function  is  a  solution  of 


(A2a) 


where  c^  is  the  dilatational  wave  speed.  If  <j>[0,  r,  (z-ct)/3^]  is  a  scalar 
wave  function,  then  <£[0,  r>  z/3^]  is  a  scalar  harmonic  function.  We  have  in 
expanded  form 


(A2b) 


We  note  that 


(A3) 
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sin  (z-ct) 


62 


Figure  Al.  Sinusoidal  Pressure  Pulse  in  a  Hollow  Cylinder 


so  that 
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+ 


1  it  + 

r  3r 


1 


2 


r 


+ 


[c 


1 


2 

-c 


-  0 
dz 


(  A4) 


But 


3 


2 

1 


We  define 


5 


1 


=  z/gj 


( A5) 


Then 


a2 4>  _i  3_i 

3r  2  r  r 


+ 


3d) 


0 


(A6) 


Analysis  of  the  vector  potential  yields  similar  results.  If  r,  (z- 

ct )/  &  ]  is  a  solution  of  the  vector  wave  equation 

2  A 

2„2;  3> 

c?  V  i|»  =  — f  (A7) 

3tZ 

then  \J>[0,  r,  a  solution  of  the  vector  harmonic  equation 


(A8) 


where  C2  is 


the  velocity  of  the  shear  wave  and 


,  2  2W  2 
(c2  "C  )/c2 
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We  assume  that 


0  <  c  <  c 


2 


(A9) 


in  the  preceding  analysis. 

The  finite  series  for  G(  0,  r,  z-ct)  can  be  derived  from  harmonic  scalar 
potentials  in  the  variables  0,  r,  z/g^  and  the  vector  harmonic  potentials  in 
the  variables  in  the  0,  r,  z/g It  is  sufficient  to  consider  scalar  and 
vector  harmonic  functions  in  the  variables  0,  r,  z,  as  the  required  potentials 
in  terras  of  0,  r,(z-ct)  can  be  obtained  by  appropriate  changes  in  the 
independent  variable  z. 

We  consider  the  scalar  form  of  Laplace's  equation  in  detail.  We  have 


(A10) 


We  assume 


<J>  =  F  sin(n0) 


(All) 


then 


32F  +  J.  3F 


2  2 

n  F  jTF 


2  2 
3r  F  3zZ 


0 


(A12) 


The  function 


F 


o 


(Aq  +  Aj  log  r)(l  +  0^) 


(A13) 
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is  an  obvious  solution  when  n  =  0*  For  n  >  0  we  have  solutions  of  the  type 


F  =  (AQr  n  +  A1rI1)(l+31z) 


(A14) 


The  general,  separable  solutions  involve  Bessel  and  trigonometric  functions. 

We  require  nonseparable  solutions  in  ascending  powers  of  z.  Both  positive  and 
negative  powers  of  r  are  required  in  order  to  satisfy  boundary  conditions  at 
the  inner  and  outer  cylindrical  surfaces.  The  required  harmonic  function  can 
be  obtained  by  expressing  spherical  harmonics  in  terms  of  cylindrical 
coordinates  and  rearranging  terms  in  ascending  powers  of  z,  References  A1 
and  A2.* 

These  potentials  can  be  obtained  in  an  elementary  manner  by  assuming 


F 


k 


G  + 
o 


4 

z 


G.  +. .  • 

4 


(A14a) 


or 


F 


k+1 


z  G, 


(A14b) 


where 


A^e.W.  Hobson,  The  Theory  of  Spherical  and  Ellipsoidal  Harmonics, 
University  Press,  Cambridge,  1931 .  See  Chapter  IV. 


A2 

J.  Dougall,  ftAn  Analytical  Theory  of  the  Equilibrium  of  an  Isotropic 
Elastic  Rod  of  Circular  Cross  Section , "  Transactions  of  the  Royal  Society 
of  Edinburg,  1913,  Vol.  XLIX,  Part  IV  (Ho.  17),  pages  895-978. 

*The  fourth  degree  harmonic  *  given  on  page  348,  Reference  12,  occurs  in  the 
torsion  problem.  Retails  will  be  described  in  a  forthcoming  BRL  report. 
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The  remaining  functions,  which  are  also  functions  of  r  alone,  can  be 
obtained  recursively  from  the  formula 


2 

q  2 

TT-  *  7  -r  -  h  Gi  “  (1+2)(i+1)  Gi+2  <A16> 

3r  r 


which  is  obtained  by  substituting  the  preceding  expressions  for  into  Eq . 
(A10),  carrying  out  the  indicated  operations,  and  equating  the  successive 
powers  of  z  to  zero. 

As  an  example,  set 


n  -  0,  k  =  4,  A  =  3 ,B  =  3 


(A17) 


We  find 


4  2  2  4 

F.  =  3r  -24r  z  +  8z 
k 


„  -4  -6  2 

+  3r  -24r  z  + 


8r~8  z4 


(  A18) 


The  potential  given  by  the  first  line  can  be  used  in  the  analysis  of  torsional 
loads  and  differs  only  by  a  constant  from  the  potential  (j>^  given  in  Reference 
12,  page  348. 

Next,  we  consider  the  contribution  of  the  Fourier  integral  to  the  total 
solution.  We  have  either 


Fg ( 9 >  r,  z) 


sin(n0 ) 


/ 


N0(r,s) 

D(s) 


sin[s(z-ct) ] 


—00 


S 


ds 


(A19a) 
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Figure  A2.  Stepped  Pressure  Pulse  in  a  Hollow  Cylinder 


or 


Fc(0,  r,  z) 


00 

cos(n0)  j 

— oo 


Na^r,S)  sins  s(z-ct)  , 
- - - -  ds 


(A19b) 


for  stresses  which  are  odd  in  z.  The  cosine  replaces  the  sine  under  the 
integral  sign  if  the  stress  is  even  in  z.  For  values  of  z  which  are  not  too 
small,  we  calculate  F  or  F  by  the  theory  of  residues,  which  is  validated  in 
cases  of  interest  by  Jordan's  lemma.  If  z>ct  we  use  a  contour  in  the  upper 
half  plane,  indented  at  the  origin  of  the  complex  s  plane.  If  z<ct  we  use  the 
corresponding  contour  in  the  lower  half  plane.  The  complex  variable  theory 
follows  suggestions  in  a  landmark  paper  by  Dougall  (A2).  Dougall  uses  an 
entirely  different  set  of  potentials  to  calculate  the  displacements  and 
strain;  however,  his  remarks  on  the  nature  and  distribution  of  the 
characteristic  roots  are  still  valid. 

It  is  obvious  the  characteristic  functions  for  the  solid  rod  do  not  have 
a  branch  point  at  the  origin  of  the  complex  s  plane,  since  only  ordinary  or 
modified  Bessel  functions  of  the  first  kind  are  involved.  These  functions  are 
analytic  in  the  entire  complex  plane  with  the  exception  of  essential 
singularities  of  exponential  type  at  infinity.  A  detailed  analysis  is 
required  for  the  hollow  cylinder,  since  functions  of  the  first  and  second  kind 
are  both  required  in  order  to  satisfy  boundary  conditions  at  the  inner  and 
outer  cylindrical  surfaces.  Bessel  functions  of  the  second  kind  have  a 
logarithmic  singularity  at  the  origin,  plus  a  finite  Laurent  series  in  the 
case  of  functions  of  integral  orders  greater  than  zero.  We  must  show  the 
characteristic  functions  do  not  have  a  branch  point  at  the  origin;  multiple 
poles  will,  in  general,  occur  due  to  the  reciprocal  powers  in  the  Laurent 
series . 

To  this  end,  we  write  D(s)  in  extended  form.  We  have  from  Eq .  (164) 
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S^iU) 

S1 ,2(a) 

si>3(a) 

Si,4(a) 

Sl,5(a) 

S1 ,6<a) 

S4,i<a) 

S4,2^a) 

S4 , 3 (a ) 

S4i4(a) 

s4j5(a) 

s4,6<a> 

D(  s 

)  =  s5 , 

S5 j2(a) 

S5 , 3 ( a ) 

S5j4(a) 

S5 , 5 (a ) 

S5,6^aJ 

Sl,2<b> 

Sl,3<b> 

Si,4(b) 

Si}5(b) 

sl,6(b) 

S4,i<b) 

s4>2(b) 

s4>3(b) 

s4,4(b) 

s4>5(b) 

s4,6<b> 

S5i<b) 

S5,2(b) 

S5 , 3 (b) 

S5>4(b) 

s5>5(b) 

s5,6<b> 

(164) 

To 

eliminate  the 

apparent 

branch  point 

from  the 

determinant  D(s)  we 

manipulat 

the 

columns  in  pairs.  Columns  1,  3,  and  5  are  derived  from  the  even-numbered 

equations  in  the 

sequence 

of  Eqs.  (80- 

115).  These  columns  therefore  contain 

only  modified  Bessel  functions  of  the 

first  kind 

,  and  are 

analytic 

functions 

of 

s.  Columns  2 

,  4 ,  and  6 

are  derived 

from  odd- 

numbered 

equations 

of  the 

above  sequence,  and  contain  only  modified  Bessel  functions  of  the  second 
kind.  The  logarithmic  term  of  the  modified  Bessel  functions  contain  the 
logarithm  of  s,  which  can  be  eliminated  algebraically The  logarithms  of  a 
and  b  which  also  occur  can  not  be  eliminated.  For  convenience  suppose  n  is 
even.  The  variable  r  in  the  following  equation  can  represent  either  a  or  b . 


K  (sr)  =  L  (s)I  (sr)  +  Ln(r)I  (sr)  +  L  (sr)  +  P  (sr) 
n  n  n  n  n  n 


(20) 


where  Ln(sr)  is  a  Laurent  series  having  a  pole  of  order  n  at  the  origin  in  the 
complex  s  plane  and  ?n(sr)  is  a  power  series  in  ascending  powers  of  (sr),  with 
zero  and  positive  powers  only.  Similarly, 


A3M.  Abramowitz  and  S.A .  Stegun,  Editors.  Handbook  of  Mathematical  Functions . 
with  Formulas ,  Graphs,  and.  Mathematical  Tables.  U.S.  Department  ofCommerce, 
Applied  Mathematic  Series  Ho.  55,  1964.  See  Eq.  (9.6.11),  page  375. 
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Vi(sr) 


K  (sr)  =  -  L  (s)I  (sr)  -  Ln(r)I  (sr) 
n+1  n  n  n 


-  p„+i<sr) 


(A21) 


Now  multiply  columns  1,  3,  and  5  by  Ln(s)  and  subtract  the  results  from 
columns  2,  4,  and  6, respectively .  We  find  the  logarithm  of  s  is  eliminated. 

If  n  is  an  odd  integer  we  add  instead  of  subtracting.  In  either  case,  we  find 
D(s)  is  an  analytic  function  of  s  except  for  poles.  A  multiple  pole  generally 
occurs  at  the  origin. 

It  is  believed  the  characteristic  roots  are  simple  except  at  the 
origin.  Asymptotic  methods  may  be  used  for  characteristic  roots  of  large 
modulus.  The  smallest  nonzero  root  can  be  examined  by  expanding  D(s)  in  a 
Laurent  series.  Extensive  calculations  have  shown  the  characteristic 
equations  arising  from  axisymmetric  loading  are  simple  for  a  large  range  of 
wall  ratios  (A4).  Moreover,  a  double  zero  leads  to  functions  of  the  type 
R(r)  x  [s(z-ct)  sin  s(z-ct)  +  cos  s(z-ct)],  which  do  not  satisfy  the  wave 
equations.  Classical  methods  of  determining  the  multiplicity  of 
characteristic  roots  are  forbiddingly  difficult  when  applied  to  the 
determinant  D(s);  hence,  detailed  analysis  will  be  reserved  for  specific  types 
of  loading  under  consideration. 


A4A.S.  Elder  and  K.L.  Zimmerman ,  "Stresses  in  a  Gun  Tube  Produced  by  Internal 
Pressure  and  Shear,"  BRL  Memorandum  Report  No •  2495,  June  1975  (AD  A012765). 
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ATTN:  DRCPM-AAH-30mm 

Dover,  NJ  07801 

2  Product  Manager 

M110E2  Weapon  System,  DARCM 
ATTN:  DRCPM-M110E2 

Rock  Island,  IL  61299 

3  Commander 

US  Army  Research  Office 

P.0.  Box  12211 

ATTN:  Technical  Director 

Engineering  Division 
Metallurgy  &  Materials 
Division 

Research  Triangle  Park 
NC  27709 

2  Commander 

US  Army  Research  Office 
ATTN:  Dr.  J.  Chandra 

Dr.  F.  Schmiedeshof f 
Research  Triangle  Park 
NC  27709 

4  Director 

US  Amy  Mechanics  and 

Materials  Research  Center 
ATTN:  Director  (3  cys) 

DRXMR-ATL 

Watertown,  MA  02172 

1  Director 

US  Amy  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

1  Commander 

Naval  Air  Systems  Command 
Washington,  DC  20360 


1  Commander 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

2  Commander 

David  W,  Taylor  Naval  Ship 

Research  &  Development  Center 
Bethesda,  MD  20084 

1  Commander 

Naval  Air  Development  Center, 
Johnsville 

Warminster,  PA  18974 

1  Commander 

Naval  Missile  Center 
Point  Mugu,  CA  93041 

2  Commander 

Naval  Surface  Weapons  Center 
Dahlgren,  VA  22448 

2  Commander 

Naval  Surface  Weapons  Center 
Silver  Spring,  MD  20910 

2  Commander 

Naval  Weapons  Center 
China  Lake,  CA  93555 

1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 

1  Superintendent 

Naval  Postgraduate  School 
ATTN:  Dir  of  Lib 

Monterey,  CA  93940 

1  Commander 

Naval  Ordnance  Station 
Indian  Head,  MD  20640 

2  AFRPL 

ATTN:  W.  Andrepont 
T.  Park 

Edwards  AFB,  CA  93523 
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AFATL 

Eglin  AFB,  FL  32542 
AFWL/SUL 

Kirtland  AFB,  NM  87117 


Director 

Lawrence  Livermore  Lab 
P.O.  Box  808 
Livermore,  CA  94550 

Director 

Los  Alamos  Scientific  Lab 

P.O.  Box  1663 

Los  Alamos,  NM  87545 


Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 

Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
Hampton,  VA  23365 

Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
ATTN:  Library 

Houston,  TX  77058 

BLM  Applied  Mechanics 
Consultants 
ATTN:  Dr.  A.  Boresi 

3310  Willett  Drive 
Laramie,  WY  82070 

CALSPAN  Corp. 

ATTN:  E.  Fisher 

P.O.  Box  400 
Buffalo,  NY  14225 


1  S&D  Dynamics,  Inc. 

ATTN:  Dr.  M.  Soifer 

755  New  York  Avenue 
Huntington,  NY  11743 

1  Southwest  Research  Institute 

ATTN:  P.  Cox 

8500  Culebra  Road 
San  Antonio,  TX  78228 

1  Stanford  University 

Stanford  Linear  Accelerator 
Center 

ALAC,  P.O.  Box  4349 
Stanford 
CA  94305 


1  State  University  of  New  York 

Department  of  Mathematics 
ATTN:  Dr.  Ram  Shrivastav 

Stony  Brook,  NY  11790 

2  University  of  Delaware 
Department  of  Mathematics 
Department  of  Mechanical  Eng; 
Newark,  DE  19711 

1  Princeton  University 

Department  of  Civil 
Engineering 

ATTN:  Dr.  A.C.  Er ingen 

Princeton,  NJ  08540 

1  University  of  Kentucky 

Department  of  Computer 
Science 

ATTN:  Prof  H.C.  Thacher,  Jr 

915  Patterson  Office  Tower 
Lexington,  KY  40506 
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2  University  of  Wisconsin- 

Mad  is  on 

Mathematics  Research  Center 
ATTN:  Dr.  John  Nohel 

Dr .  Richard  Meyer 
610  Walnut  Street 
Madison,  WI  53706 

1  Virginia  Commonwealth  Univ. 

Department  of  Math.  Sciences 

901  W.  Franklin  St. 

Richmond,  VA  23284 


1  University  of  Colorado 

Dept,  of  Aerospace  Eng.  Sciences 
ATTN:  Prof.  G.  Inger 
Boulder,  CO  80304 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP ,  H. 

Dir,  USAMTD 

ATTN:  H.  King 

P.  Paules 

Cdr ,  USATECOM 

ATTN:  DRSTE-TO-F 

Dir,  USACSL, 

Bldg.  E3516 ,  EA 
ATTN:  DRDAR-CL 

DRDAR-CLB 
DRDAR-CLB-PA 
DRDAR-CLD 
DRDAR-CLN 
DRDAR-CLY 
DRDAR-CLJ-L 


Cohen 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


